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ABSTRACT 

Context. Rapidly rotating, chemically homogeneously evolving massive stars are considered to be progenitors of long gamma-ray 
bursts. 

Aims. We present numerical simulations of the evolution of the circumstellar medium around a rapidly rotating 20 M© star at a metal- 
licity of Z = 0.001. Its rotation is fast enough to produce quasi-chemically homogeneous evolution. While conventionally, a star of 
20 M© would not evolve into a Wolf-Rayet stage, the considered model evolves from the main sequence directly to the helium main 
sequence. 

Methods. We use the time-dependent wind parameters, such as mass loss rate, wind velocity and rotation-induced wind anisotropy 
from the evolution model as input for a 2D hydrodynamical simulation. 

Results. While the outer edge of the pressure-driven circumstellar bubble is spherical, the circumstellar medium close to the star 
shows strong non- spherical features during and after the periods of near-critical rotation. 

Conclusions. We conclude that the circumstellar medium around rapidly rotating massive stars differs considerably from the sur- 
rounding material of non-rotating stars of similar mass. Multiple blue-shifted high velocity absorption components in gamma-ray 
burst afterglow spectra are predicted. As a consequence of near critical rotation and short stellar evolution time scales during the last 
few thousand years of the star's life, we find a strong deviation of the circumstellar density profile in the polar direction from the \/R^ 
density profile normally associated with stellar winds close to the star. 
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• 1. Introduction 

Within the collapsar picture (Woosley ' 1993K long gamma-ray 
, bursts (ORB) are thought to be produced at the end of the evolu- 
' tion of massive stars, if their rapidly rotating iron cores collapse 
to form black holes. The circumstellar medium (CSM) around 
ORB progenitors is interesting because it determines the ORB 
afterglow evolution (Chevalier & Li l200QI Panaitesku & Kumar 
I2QQ21 12QQ2I Chevalier et al. 20041 Ramirez-Ruiz et al. 2005. 
Nakar & Granot 2006 ), and because it may produce distinct ab- 
sorption features in GRB afterglow spectra (Mirabal et al. 120031 
Schaefer et al. 2003, Fiore et al. 2004, Starling et al . [2005]). In 
previous papers, van Marie et al. (i2005bl I2QQ6I I2QQ71) have com- 
puted the CSM evolution of various types of supernova and GRB 
progenitors without including the effect of rotation of the central 
star. 

It has been recently suggested that rapidly rotating, chemi- 
cally homogeneous massive stars evolve into long GRBs (Yoon 
& Langer '200 3 Wo osley & Heger '2005^; Yoon et al. 2006; 
Cantiello et al. 12007 1. This type of evolution is thought to oc- 
cur at low metallicity, where stellar winds are weaker and an- 
gular momentum loss by winds (Langer . 1998 J is therefore lim- 

Send offprint requests to: A. J. van Marie 



ited. Since models of chemically homogeneously evolving stars 
avoid a redward evolution in the HR diagram, but rather evolve 
from the main sequence directly to the helium main sequence, 
the structure of their circumstellar medium is expected to be dif- 
ferent from that of Galactic Wolf-Rayet stars, which are thought 
to be the descendants of red supergiants or Luminous Blue 
Variables (Garcia-Segura et al.' 1996al[r996bb . Additionally, the 
chemically homogeneous stellar models remain very rapid rota- 
tors throughout their lives, and encounter various stages of near- 
critical rotation (Yoon & Langer 2005). Wind anisotropics may 
therefore play an important role in shaping the CSM of these 
stars. 

Here, we undertake an effort to understand the evolution of 
the circumstellar medium around such stars. We perform de- 
tailed numerical simulations of its time evolution, were we use a 
detailed stellar evolution model as input. We employ a prescrip- 
tion of the anisotropy of the mass outflow and wind velocity as 
a function of the rotation of the central star, and adopt this, to- 
gether with the time dependent global wind properties and ioniz- 
ing photon fluxes, as the inner boundary condition for our CSM 
evolution model. As result, we obtain the density, velocity and 
temperature structure of the CSM of our stellar model, from the 
main sequence to the pre-SN stage. 
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Fig. 1. Evolution the rapidly rotating 20 M© model at Z = 
0.001 which is used as input for the CSM simulations in the 
Hertzsprung-Russel diagram. The model is similar to the ones 
used in Yoon et al. (2006). Its initial rotational velocity is 
414 km/s, which is about 60% of the critical velocity. The 
track starts at the zero age main sequence and ends at the pre- 
supernova stage. 
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Fig. 2. Latitudinally averaged mass loss rate, and D (= Vrot/ Vcrit) 
for our 20 M© sequence, as a function of time. The horizontal 
axis shows the logarithm of the time before core collapse (tec)- 
The gray area highlights the time interval which is explored in 
Figs.[5]to[7]below. 



In Sect. |2] we describe the evolution of the properties of the 
chemically-homogeneously evolving stellar model which is used 
as input for our hydrodynamic calculations. In Sect. [3] we discuss 
the numerical method used to simulate the hydrodynamical evo- 
lution of the CSM. The results of our simulation are presented in 
Sect. 14. 1[ In Sect. [5l we show the absorption profiles that would 
result from the morphology of the circumstellar medium when 
the star explodes. Section [5] deals with the density profile of the 
circumstellar medium, which is important for the lightcurve of a 
GRB afterglow. In Sect. [7] we discuss our results. 

2. The stellar evolution model 

If a star rotates very rapidly, chemical mixing may occur on such 
a short timescale that no significant chemical gradient can be 
established (Maeder 1987). Since new material is mixed into 
the core at all times, core hydrogen burning in these stars will 
continue until all hydrogen in the star is exhausted. Due to this 
mixing, the star becomes a Wolf-Rayet star during the core- 
hydrogen burning phase and never swells up to become a red su- 
pergiant. At core hydrogen and core helium exhaustion, the star 
contracts to reach the necessary temperature for the next burn- 
ing stage. As the star shrinks, its rotation rate increases due to 
angular momentum conservation. Once the star reaches critical 
rotation, it can not spin-up any further. Instead, it loses mass at 
a highly increased rate. With the mass it also loses angular mo- 
mentum, until the spin velocity drops below the critical rotation 
velocity. 

At solar metallicity, massive stars have such a high mass loss 
rate that they lose large amounts of angular momentum (Langer 
119980 . As the star is spun down, the mixing time increases such 
that a chemical gradient can be established between the core 
and the envelope. However, at low metallicity, the mass loss 
rate even of Wolf-Rayet stars is considerably lower (Nugis & 
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Fig. 3. Number of high energy photons, and terminal wind ve- 
locity of our 20 M© sequence star, as a function of time. The 
horizontal axis shows the logarithm of the time before core col- 
lapse (tec)- The gray area highlights the time interval explored in 
Figs.[5]to|7]below. 



Lamers|2000l Crowther et al. [2002 : Vink & de Koter |2005]), 
and it is possible that the star remains quasi-chemically homoge- 
neous throughout core h ydrogen burning (Yoon & Langer .2005t 
Woosley & Heger[2005]). 

As an input for our hydrodynamic CSM study, we use the 
results of a stellar evolution calculation for a 20 M© star with 
a metalHcity of Z = 0.001 (cf. Yoon et al. l2006l) . The initial 
equatorial rotational velocity of the model is 414 km/s. The input 
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Fig. 4. The result of a ID hydrodynamic simulation of the CSM 
around our 20 M© star, just before the star reaches critical ro- 
tation for the first time (4.644 x lO^yr before core collapse; 
cf. Fig. 2). The solid line shows the mass density, the dashed 
line the internal energy density. Moving outward from the star, 
we first encounter the free streaming wind, then the shocked 
wind material (6 < R < 30 pc), the photo-ionized ISM 
(30 < R < 43 pc), the shell of shocked ISM (7? ^ 45 pc) 
and finally, the unperturbed ISM. The dotted line marks the 
Stromgren radius. We use this point in time to start the 2D sim- 
ulation. 



physics in the stellar evolution calculation are as in Yoon et al. 
(I2006K but without stellar wind enhancement due to increased 
CNO surface abundances. Fig. [T] shows the evolutionary track 
of our star in the Hertzsprung-Russel diagram; as a result of the 
quasi-chemically homogeneous evolution, the star evolves from 
the main sequence directly into the Wolf-Rayet regime. 

The mass loss history of chemically homogeneously evolv- 
ing stars is diff'erent from that of non-rotating or slowly rotat- 
ing stars. The latter have a low mass loss rate during their main 
sequence, and a high mass loss rate during the red supergiant 
stage. The mass loss history of our homogeneous star is shown in 
Fig. [21 Mass loss is low, except during two short phases. During 
these phases the star reaches critical rotation (dashed line in Fig. 
|2]). This happens for the first time during the overall contrac- 
tion of the star following core hydrogen exhaustion. During this 
stage, the star loses about 2 M© in an equatorial outflow, in order 
to avoid overcritical rotation. Upon core helium exhaustion, the 
star shrinks again toward carbon ignition, and approaches crit- 
ical rotation again. However, this time the star reaches the end 
of its evolution and explodes as a supernova before its rotational 
velocity can drop again. The amount of mass lost in the second 
phase of near-critical rotation is 0.05 M©. 

The terminal wind velocity used for our hydrodynamical 
calculations does not follow directly from the stellar evolution 
model. We assume here that the wind is radiation driven, which 
leads us to the following formula 



e.g. Cherchnefl' & Tiellens (119941) , where Fe = L/Le, the ratio 
between stellar luminosity and Eddington luminosity and /3 is a. 
parameter related to the surface temperature of the star (Eldridge 
et al. 120051) . The resulting wind velocity is shown in Fig. [3] Due 
to the compactness of the stellar model throughout its evolution, 
the wind speed is at all times of the order of a few thousand 
kilometers per second. However, pronounced minima occur at 
the times of near critical rotation. 

The values for the mass loss rate and wind velocity shown 
in Figs. [2] and [3] are spatial averages. Rapidly rotating stars have 
winds that depend strongly on the latitude. As long as the rota- 
tional velocity is well below critical, this efl'ect is not important. 
However, as the rotational velocity increases, so does the as- 
phericity of the wind. We use the equations found by Bjorkman 
& Cassinelli ( 1993, from here BC) for winds from rapidly ro- 
tating stars (see Sect. 3). This model assumes that for a rapidly 
rotating star, the wind will be denser at the equator than at the 
poles, while the wind velocity will be higher at the poles than at 
the equator (cf. Sects. [3l and 14. II) . 

Since we include the efl'ects of photo-ionization on the CSM, 
we have to find the number of high energy photons emitted by 
the star. We approximate this by treating the star as a black body 
radiator. This means that the number of photons emitted per sur- 
face area with a certain wavelength only depends on the surface 
temperature of the star, according to Planck's formula. We sim- 
ply integrate the total number of photons above the hydrogen 
ionization threshold (13.60 eV) and correct this number for the 
diminishing photo ionization cross section at higher energies as 
given by Osterbrock (119741 . The resulting number of ionizing 
photons radiated by the star are shown in Fig.O 

3. Numerical method 



We use the ZEUS 3D code (Stone & Norman[T992l Clark[l996| 
for our simulations. Our numerical method is the same as used 
by Garcia-Segura et al. (1996a & 1996b) and in our previous 
papers (van Marie et al. 2005a, 2005b and 2007). This means 
that we start the simulation in one dimension and map the result 
onto a two dimensional grid when it becomes necessary. Even 
in the case of a rotating star it is possible to do this. Although 
the wind is aspherical, the main sequence shell is driven by the 
thermal pressure of the shocked wind material (Weaver et al. 
119771) . Therefore, the outer shell of the main sequence bubble 
will still be spherical, rather than ellipsoid. The wind termina- 
tion shock of course is not spherically symmetric, since the ram 
pressure of the wind depends on the polar angle. However, this 
is no problem, if we start the 2D simulations in time for the gas 
to compensate for this change. 

In the 2D simulations we use the same equations for the an- 
gle dependence of the wind parameters at large distances from 
the star as were used by Langer et al. (1999) and Chi{a et al. 
( 2007) . These are based on the model in BC. So the wind veloc- 
ity is: 



2A/3GM(1-Te) 



R 



while the wind density is: 

(a/2)BM(l-^sinO)^ 



P(0) 



(2) 



(3) 



Voo = 



2ySGM 



R 



where we set the parameters defined in BC to y = 0.35 and 
(1) ^ = - 0.43. The eff'ect of the rotational velocity is included in 
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the parameter Q = Vrot/ Vcrit- Finally, the parameter a is defined 
by: 



a = 



9 / . Q sm ^ , , . 

cos + cot ^ 1 + y : — sm 

' 1 - 12 sm 



-1 



(4) 



with 0' = Q sin6/Vcrit/[2 ^/2 Yoo(0)]. Note that for Q > 0.995 
these equations do not have a valid solution, since 0' becomes 
larger than n/2, which would start to reverse the process of con- 
centrating the mass loss toward the equator. Therefore, we set 
Q < 0.99 in our simulations. The correction factors A and B 
serve to make sure that the average stellar mass loss and wind 
velocity remain the same as in the ID approximation, by averag- 
ing the latitude dependent parts of Eqs. [2] and [3] over the surface 
area of the star and normalizing this to 1 . They are calculated by 
integrating Eqs. [2l and [3] over the surface of the star and normal- 
izing the result to the average values shown in Fig. [2] and O For 
a diff'erent approach to the problem of winds from rotating stars 
and the eff'ect on the CSM see Eldridge (l20Q7l) . 

In our simulations we include the eff'ect of photo-ionization 
in the same way as van Marie et al.( 2005b, 2007 ). This means 
that we calculate the Stromgren radius along each radial grid 
line. All matter within this radius is presumed to be fully ion- 
ized, while all matter outside the Stromgren radius is consid- 
ered neutral. In a completely realistic model, the photon count 
would depend on the latitude, since the star is no longer spher- 
ical at high rotational velocity. However, since our treatment of 
the photo-ionization is only an approximation, we treat the pho- 
ton field emitted by the star as if it were spherically symmetric. A 
more sophisticated way of treating the photo-ionization by mas- 
sive stars was presented by Arthur & Hoare ( 2006) and Mellema 
et al. (l2QQ6b . 

The star has two periods in which the wind is strongly as- 
pherical. The first lasts about 10000 years and the mass shed 
during this phase will have enough time to spread throughout 
the entie circumstellar bubble. However, the second period is 
very short and is immediately followed by the supernova explo- 
sion. Therefore, its eff'ects are only visible close the star. This 
means that we have to perform a special simulation of the CSM 
close to the star for the last few thousand years of its evolution. 
Fortunately, due to the short period followed immediately by 
core collapse, the eff'ects of this second period can only be felt 
close to the star. Therefore, this new simulation only has to cover 
the free- streaming wind region. The hot bubble of shocked wind 
material, which is typical of wind blown bubbles (Weaver et al. 
11977 ') does not have to be taken into account. This makes the 
calculation easier since a simple outflow boundary can be used 
at maximum radius rather than the more complicated construc- 
tion used by Garcia-Segura et al ( 1996a, 1996b), for a situation 
where the outer boundary was inside the hot bubble. 

4. Results 

4.1. Core hydrogen burning 

For most of the core hydrogen burning, the star rotates at about 
= 0.6 (Fig. [2]), which does not produce significant aspherici- 
ties. Moreover, even though the wind termination shock may be 
slightly aspherical, the main sequence bubble is pressure driven 
and will thus be very close to a spherical shape anyway. For 
this reason, we modeled the evolution of the CSM of our 20 M© 
star in ID, up to shortly before the first sharp rise in Q (Fig. [2]). 
For this we use a spherically symmetric grid, with 2000 ra- 
dial grid points. The Interstellar medium (ISM) density is set 




Fig. 8. Density distribution of the circumstellar medium 1 .47 x 
lO^yr before the star explodes as a supernova (cf.. Figs. [5] to [7]). 
The eff'ects of the first critical rotation phase have long since dis- 
appeared. The number of ionizing photons has decreased, while 
the ram pressure of the wind is larger than before the first period 
of critical rotation. Together, these factors cause the hot bubble 
of shocked wind material to expand at the expense of the sur- 
rounding H II region. The eff'ect of the second period of critical 
rotation cannot yet be seen, since this becomes visible only in 
the last few hundred years before the supernova (See Sect. 14.31) . 



to 10"^^-^g/cm^. The result of this simulation is shown in Fig. 
m This result closely resembles that of a non-rotating 25 M© 
star (van Marie et al. 2004). Moving outward from the star, we 
first encounter the free- streaming wind, then the shocked wind 
material, then the HII region of photo-ionized ISM, the shell 
of shocked ISM material and finally the unperturbed ISM. As 
in van Marie et al. (|2004| ), the hot shocked wind material and 
photo-ionized ISM have reached pressure equilibrium. This is 
not always the case, since a strong wind can sweep up the entire 
H II region, as was shown in Freyer et al. (120031) and van Marie 
et al. (2006). The star reaches critical rotation for the first time 
after about 1.2868x10^ years (Fig. [2]). Therefore we start out 2D 
simulation about 30 000 years before this time. 

4.2. Core hydrogen exhaustion through heiium burning 

To switch from ID to 2D, we have mapped the result of the ID 
simulation onto a 2D grid, which has 1000 radial and 200 angu- 
lar grid points for an angular size of tt/ 2 (Fig. [5]). 

As the star reaches critical rotation, the mass loss rate in- 
creases by nearly three orders of magnitude (from 10"^ to 
10"^ M© yr"^), while the average wind velocity decreases by 
about 600 kms"^ This eff'ect is largest for near-equatorial lat- 
itudes. Since the ram pressure of the wind varies considerably 
with latitude, the wind termination shock becomes very aspher- 
ical. Due to the high mass loss rate, a new shell is formed at the 
wind termination shock. Inside the wind termination shock, the 
free- streaming wind becomes strongly anisotropic, with most of 
the mass flow concentrated close to the equatorial plane (Fig. [5]). 
since the 



van Marie et aL: Circumstellar medium around a chemically homogeneous star 



5 



The period of critical rotation only lasts for about 
10000 years. Afterwards, the mass loss rate decreases and the 
wind becomes almost spherically symmetric once again. The 
Wolf-Rayet wind sweeps up the disk, forming another shell, 
which moves outward (Fig. [6]). As a result, the aspherical shell 
at the original wind termination shock is broken up and the rem- 
nants of both shells are pushed out into the hot bubble, where 
they will eventually dissipate. The point where the wind breaks 
through the shell depends on the exact wind parameters dur- 
ing and after the period of critical rotation. Effectively, the disk 
around the equator deflects part of the wind upward, which in- 
creases the ram pressure at higher latitude. At the same time, the 
shell is thinner at higher latitude, since the mass loss rate was 
lower in that direction during the period of critical rotation. By 
the time the star explodes as a supernova (t ~ 1.331 x 10^ yr), 
the remnants of the shell have completely disappeared (Fig. [5]). 

In the phase between the two periods of critical rotation (the 
core helium burning phase of the star), the corrected number of 
ionizing photons is smaller than during the main sequence phase, 
while the ram pressure of the wind increases. Due to the lower 
number of photons, the Stromgren radius recedes, leaving the 
original main sequence shell (R ^ 45 pc) unsupported. This 
shell starts to dissipate, while a new shell is formed at the transi- 
tion between shocked wind material and the former H II region. 
This shell, driven by the thermal pressure of the shocked wind 
material, starts to expand into the former H II region. This expan- 
sion goes quite fast, since the density in the former H II region 
is lower than in the unperturbed ISM. A similar efl'ect is found 
for the CSM around non-rotating massive stars during the Wolf- 
Rayet phase (van Marie et al. 2005b). 

If the left hand image in Fig. [6] were to be rotated around the 
equator, its shape resembles the ring structure observed in SN 
1987 A. Since it is a temporary feature which quickly disappears, 
it may indicate that the rather peculiar shape of SN 1987 A is the 
result of a period of rapid rotation of the progenitor star some 
time in the recent past. 

4.3. Beyond core helium burning 

The stellar wind evolution during the final ^ 2000 yr of the star's 
life as shown in Figs. [2l and [3] produces dramatic changes in the 
CSM in the inner few parsecs around the central star. In Fig. 
[SI the final frame of our first 2D-hydro run, some non-spherical 
features can be seen close to the star, while the rest of the CSM 
remains spherically symmetric. The eff'ects of the first phase of 
near-critical rotation have disappeared. The further evolution is 
followed by a second 2D simulation with higher time and spacial 
resolution, and restricted to the inner 5 pc of the CSM. Due to 
the short remaining time to core collapse, the outer parts of the 
bubble will remain largely unchanged. 

Just before the supernova explosion, the star reaches critical 
rotation for the second time (Fig.©. This period lasts only about 
1000 yr, so no shell is formed at the wind termination shock, and 
the eff'ects of the increased rotational velocity are not felt at large 
radius (see also Sects. [3] and 14. lb . However, close to the star, the 
wind is aspherical once again. In Figs. [9] and \TU\ we show the 
mass loss rate, fl, the high energy photon count, and the wind 
velocity during the last stage of the evolution of the star. 

The evolution of the CSM close to the star during the last 
500 years is shown in Fig.[TT] Eff'ectively, it repeats the evolution 
of the first stage of critical rotation. At first, the wind velocity de- 
creases (Fig. [9]). This creates a rarefaction wave, which moves 
outward in the free- streaming wind. Then the wind becomes 
strongly aspherical and the mass loss rate increases (Figs. [9] and 
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Fig. 10. Number of ionizing photons, and D (= Vrot/ Vcrit), as a 
function of time, during the same phase as shown in Fig.[9l We 
use the average photon flux for all angles. 

fTOl) . At the same time, the wind velocity increases again. A shell 
starts to form at the inner edge of this rarefaction wave, due to 
the increased mass loss rate and the fact that the new, faster wind 
starts to sweep up the material in the rarefaction wave. However, 
before the shell has an opportunity to fully form, the evolution 
of the star is over and it explodes as a supernova, leaving Fig.fTTI 
as the final morphology of the CSM close to the star, into which 
the GRB and supernova will expand. 

The 1 /R^ density distribution of the free streaming wind is 
compromised, not only by the rarefaction wave, but also by small 
changes in mass loss rate and wind velocity during the final years 
of the stellar evolution. 



5. Velocity-resolved column depth 

In the afterglow of GRB 021004 a number of blue-shifted ab- 
sorption lines were observed, indicating the presence of matter 
moving away from the progenitor star at discrete velocities, van 
Marie et al. (12005b ) showed that these absorption features might 
be explained by the hydrodynamic interactions in the circum- 
stellar medium, as proposed by Schaefer et al. (2003 ), Mirabal 
et al. (2003 ), Fiore et al. (2004) and Starling et al. (2005 ). Since 
the rapidly rotating star has a diff'erent evolution from that of a 
non-rotating star, we expect to see a diff'erent absorption pattern 
here. 

In order to calculate the blue-shifted absorption pattern we 
use the same method as in van Marie et al. (2 005b[|2007b . This 
gives us the column density per 1 km/s interval as a function of 
radial velocity, which corresponds to the absorption as a func- 
tion of blue-shift. We correct for the Maxwell-Boltzmann distri- 
bution of particle velocities according to the temperature of the 
gas. It is clear from Fig. [8] that the first phase of critical rotation 
will have no observable consequences at the time of the star's 
death. Therefore, we concentrate on the CSM interactions due 
to the second period of critical rotation, described in Sect. 14. 3[ 

While van Marie et al. ( 2005b, 2007 ) produced the absorp- 
tion pattern as a function of time averaged over all angles, this 
is not useful here due to the strong asphericities in the CSM. 
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For the GRB context, we are primarily interested in the poten- 
tial absorption pattern in a GRB afterglow, which means that we 
focus on the material that lies along the polar axis of the star. 
Therefore, we average the absorption pattern found along the 
first 20 radial grid lines, which corresponds to a cone with an 
opening angle of 18° around the polar axis (Fig.[T2l). 

Figure [T3] shows the column density as a function of radial 
velocity for the final 1500 yr of the evolution of the star. In order 
to consider also the material beyond a radius of 5 pc, we have 
combined the final outputs of both 2D simulations as follows. 
For the inner 5 pc, we use the final output shown in Fig.[TTl For 
the rest of the circumstellar bubble, we use the output shown in 
Fig.IH 

A strong component at zero velocity is always visible. This 
originates from the ISM, the swept up shell of ISM material and 
the hot bubble, where the material in the hot bubble has such 
a low density that it hardly contributes at all. Moreover, it has 
a high temperature, which smears its contribution out over an 
extremely large velocity interval. The second feature occurs at 
high velocity and changes rapidly over time. 

At the start of Fig. [13] we only see a single high velocity 
component, since the stellar wind has a nearly constant velocity 
for some time. The velocity then evolves according to Fig. \TU\ 
First it drops almost instantaneously at ^ ~ 1.330595 x 10^ yr. 
After this, the velocity increases again. This is especially true 
at the pole, since the rotational velocity increases also, and the 
stellar wind velocity for our rotating star is highest at the pole 
(Eq. O. Therefore, at the time of core collapse, we see three 
high velocity components at distinct velocities: 

1. The original wind velocity, which still exists in the free- 
streaming wind further away from the star. 

2. The lowest velocity reached by the wind during this episode, 
which exists in the rarefaction wave that is moving outward. 

3. The latest wind velocity, which exists in the wind closest to 
the star. 

The first two components become weaker over time, since the 
material with the original wind velocity enters the hot bubble 
where it is diluted and the rarefaction wave is being swept up by 
the faster wind that succeeds it. 

According to this model, an observer might see a double or 
triple absorption feature at high velocity. A double absorption 
feature was in fact observed for GRB 021004 (Fiore et al. ,20041 
Starling et al. 2005 ). Absorption features at intermediate veloci- 
ties (100..600 km/s), which were also observed in GRB 021004, 
are not produced in our model. For those to occur, the wind ve- 
locity in the phase of near-critical rotation would have to drop to 
well below 1000 km/s, and an Homunculus-type nebula might 
have been produced (Langer et al. 1999 ). Whether such low wind 
velocities can be produced in the final stage of chemically ho- 
mogeneously evolving stars of diff'erent mass and/or metallicity 
will be investigated in the near future. The ionization state of the 
circumstellar medium is not taken into account in our calcula- 
tions. If the gas were ionized to a high degree, spectral features 
like C IV and SilV would not be visible at all (Prochaska et al. 
12006, Chen et al. 2007 ). Therefore the spectral hnes observed 
in GRB 021004 would be produced far away from the pro- 
genitor star, indicating an uncommonly large free-streaming 
wind region, or be unrelated to the star. Alternative expla- 
nations for this contradiction may be the presence of dust in 
the stellar wind, which influences the ionization state of the 
gas, or a structured GRB beam (Starling et al. 20051, which 
would produce an afterglow over a much broader front than 
the initial GRB 




Fig. 12. To produce Fig. [131 we average the the "absorption spec- 
tra" over each radial grid line in a 1 8° cone around the polar axis. 
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Fig. 13. Column density per radial velocity interval of 1 km/s, as 
a function of radial velocity and time, for the final 1500 yr of the 
evolution of the star. The plot ends at the time of core collapse. 
The horizontal axis shows the time prior to core collapse (tec)- 
The changes in wind velocity during the last thousand years of 
the evolution produce a complex spectrum of high velocity ab- 
sorption components. 

6. The circumstellar density profile 

The density of the circumstellar medium directly influences the 
shape of the GRB afterglow light curve. Numerical and analyt- 
ical models show that many GRBs seem to occur in a medium 
that does not conform to the 1 /R^ density profile that one would 
expect for a free streaming wind of a massive star (Chevalier 
& Li 2000; Panaitescu & Kumar [20021 [20021 and ChevaHer et 
al. 2004). With a = dlogp/dlogr, small values of a can be 
explained by assuming that the GRB jet runs into the shocked 
wind, rather than into the freely expanding wind where the den- 
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sity is close to constant (a = 0). However, this requires the wind 
termination shock to be very close to the star (Wijers lTOQli van 
Marie etal. 2006, 2007). 

The location of the wind termination shock depends on the 
balance between the ram pressure of the wind and the thermal 
pressure of the hot wind bubble that contains the shocked wind. 
The bubble pressure depends primarily on the mechanical lu- 
minosity of the stellar wind during all the preceding stages of 
the stellar evolution. Therefore in order to bring the wind ter- 
mination shock close to the star, the ram pressure of the wind 
during the final stages must be relative low, compared to the ram 
pressure of the wind during earlier phases. For a normal Wolf- 
Rayet type star, this is difficult to achieve, since Wolf-Rayet star 
winds tend to have both high velocity and a high mass loss rate. 
However, the chemically quasi-homogeneous model described 
in this paper has a very low metallicity and a completely dif- 
ferent evolution. As a result of the lower metallicity, the Wolf- 
Rayet wind is not as strong as for solar metallicity stars. The 
thermal pressure in the shocked wind bubble is not very diff'er- 
ent form the high metallicity case, since it depends on the me- 
chanical luminosity of the wind to the power 0.4 (Weaver et al. 
119771 ). Therefore, the radius where those two pressures are in 
balance will be closer to the central star for low metallicity stars. 
A second factor is the lack of a red supergiant phase. Initially, 
a red supergiant wind constrains the wind termination shock of 
the Wolf-Rayet wind, which has to plow through the dense ma- 
terial of its predecessor (Garcia-Segura et al. |1996bl van Marie 
et al. I2005bl) . However, once the Wolf-Rayet wind driven shell 
reaches the hot bubble created by the main sequence wind, a 
new equilibrium is reached between the thermal pressure in this 
bubble and the Wolf-Rayet wind ram pressure. Since the red 
supergiant wind has a low mechanical luminosity (due to low 
wind velocity) it contributes very little to the thermal pressure 
of the hot bubble while still giving it time to cool and expand, 
causing the thermal pressure to drop quickly, which places the 
wind termination shock further away from the star. This means 
that, unless the GRB occurs before the Wolf-Rayet wind reaches 
the hot main sequence bubble, which typically happens within 
about 25 000 years (Garcia-Segura et al. 1996b, van Marie et al. 
I2005bl) , a red supergiant phase will make it more difficult to cre- 
ate a constant density medium close to the star. Since chemically 
near-homogeneous stars do not have a red supergiant phase, their 
wind termination shocks will tend to be closer. In Figs. [7] and [8] 
we can see that the wind termination shock lies approx. 5-10 pc 
from the progenitor star. This is about 30-50% of the radius of 
the wind termination shock for a 'normal' Wolf-Ratet star in a 
similar environment (van Marie et al. 120061) . In itself, this is not 
enough to explain the presence of a constant density medium in 
the GRB afterglow but it can be combined with external factors, 
such as stellar motion and high density interstellar medium to 
bring the wind termination shock close enough to the star (van 
Marie et al.|2006]). 

As evident from Fig. [8l the final density profile in the inner 
few parsec of the CSM of our 20 M© star is not well described by 
ana = -2 law. In Fig. [14] we show the final density profile along 
the polar axis produced by our simulations. This density profile 
is obtained as the average over a cone with an opening angle of 
eighteen degrees around the polar axis, equal to the angle used 
in Figs. [12] and [13] While the density follows a = -2 up to 
about R = 0.2 pc, the density slope steepens considerably, up 
to a = -3 beyond that radius, up to the steep density drop at 
7? ^ 3 pc caused by the rarefaction wave (see Fig.[TQb. 

While we can not generalize from the single case that we 
investigate here, our results show that strong deviations from a 
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Fig. 14. Density profile average over a cone with opening angle 
18 degrees around the polar axis, at the time of core collapse. 
The density profile does not conform to a distribution, as 
expected from a constant wind. Density profiles with the density 
varying as R~^, R~^, and R~^ are shown for comparison. This 
plot also shows the rarefaction wave caused by the sudden drop 
in wind velocity just prior to the period of critical rotation (See 
Fig. [TO]). 

density slope of a = -2 may be expected even in the regions 
very close to the GRB progenitor, well inside the classical wind 
termination shock. The reason is the short evolutionary time 
scale of the star itself, after core helium exhaustion, and the fact 
that near-critical rotation is reached in this phase, which greatly 
amplifies the sensitivity of the wind properties to the changes 
in the global stellar properties. While in our example case we 
found a decrease of the density slope in the inner few parsec as 
consequence of an increased wind speed in the last ~ 1000 yr, 
cases where the rotation of the stellar surface drops well below 
Q ^ 1 during the last few thousand years of the star's life are 
also found in the GRB progenitor evolution grid of Yoon et al. 
( 2006), where a values closer to zero may be expected. 

7. Conclusions 

The time evolution of the mass loss rate and wind velocity of 
a rapidly rotating, chemically homogeneous massive star diff'ers 
radically from the evolution of a non-rotating star. This diff"er- 
ence is reflected in the evolution of the CSM. As was shown 
in Sects. 14.11 and 14.31 our rapidly rotating 20 M© star pro- 
duces a large bubble, which resembles that of non-rotating stars. 
However, during the two periods of near-critical rotation, large- 
scale deviations from spherical symmetry are produced. While 
the outer edge of the bubble remains spherical, closer to the star 
the asphericity of the wind creates unique morphologies. As the 
second phase of near-critical rotation occurs immediately be- 
fore the supernova explosion, the supernova itself, and a pos- 
sible GRB, will expand into such an aspherical medium. This 
will clearly influence the supernova remnant evolution, and the 
GRB afterglow light curve. 

For a GRB, the efl'ect of a non-spherical CSM can be felt in 
two ways. First, the relativistic jet will not expand in a normal. 
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free streaming wind, but in a medium that has a far more com- 
plex structure (Sect. 14.3b . In particular, the density slope in the 
innermost few parsec, well inside the wind termination shock, 
may differ significantly from 7?"^, and may also show significant 
changes as function of radius, bf As the changes of the stellar 
wind during the last few 10000 yr depend on the internal struc- 
ture of the star at this stage, it requires the consideration of a 
large grid of stellar evolution models to find out what the typi- 
cal final density profiles would be, and how large their expected 
variation is 

Second, the medium in front of the jet does not have a con- 
stant velocity. Assuming that not all of this material is ionized 
by the GRB, it may be observable in the form of blue- shifted ab- 
sorption lines (Sect. O. Our investigation shows that the wind- 
wind interactions during the final stages of the stellar evolution 
produce multiple absorption lines at high velocity. Such lines 
were observed in the afterglow of GRB 021004. While it is de- 
bated as to whether these absorption lines are in fact circumstel- 
lar in origin (Prochaska et al. 120061 Chen et al. 12007b , the pos- 
sibility of producing them in this way is shown here for the first 
time. Furthermore, such lines might also be observed in early- 
time spectra of nearby Type Ib/c supernovae at low metallicity. 
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Fig. 5. Density of the CSM at the beginning of the 2D simulation (4.565 x lO^yr before core collapse). The morphology of the 
CSM on the left is still similar to that shown in Fig.lH However, the wind termination shock has become slightly aspherical due to 
the angular dependence of the ram pressure of the wind. The Stromgren radius lies outside these figures on the inside of the main 
sequence shell (R ^ 45 pc). In the second frame (4.327 x 10^ yr before core collaspe), the star has been at critical rotation for 
some time. Close to the star the mass of the wind is concentrated around the equator. However, the wind termination shock reaches 
its maximum radius at the poles, since the wind velocity reaches its maximum in that direction. The increased mass loss rate has 
created a shell at the inner boundary of the hot bubble. 
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Fig. 6. Similar to Fig. [3 after the period of critical rotation has passed. On the left (4.169 x lO^yr before core collaspe), the star 
no longer rotates at critical velocity and the wind is closer to spherical symmetry. At this time the wind, which has become much 
faster, is in the process of sweeping up the surrounding matter. This is easiest toward the poles, since the disk around the equator 
contains most of the mass and is more difficult to move. Note the double shell at the outer edge of the free- streaming wind. The 
outer shell (R ^ 15 pc at the equator) was caused by increased mass loss and lower wind velocity during the period of critical 
rotation. The inner shell (R ^ 8 pc at the equator) is the result of the later, faster wind sweeping up the disk. Eventually, the wind 
breaks out of the shell, which was formed at the wind termination shock. The latitude at which this happens depends on the mass 
and velocity distribution, both during and after the period of critical rotation. The disk-like structure around the equator deflects the 
wind of the next phase upward, causing increased ram pressure at higher latitude. Eventually the disk itself is swept up by the wind 
as can be seen on the right (4.089 x lO^yr before core collapse). The pressure equilibrium between the ram pressure of the wind 
and the thermal pressure of the shocked wind bubble changes rapidly during this period. This causes the wind termination shock to 
move quickly. 
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Fig. 7. Similar to Figs. [5] and [6l The shell has been completely fragmented and a new wind termination shock is being established 
between the free streaming wind and the shocked wind material (simeq 3 pc in the figure on the left, 3.852 x lO^yr before core 
collapse). The shocked wind material is highly turbulent. In the right hand frame (3.614 x lO^yr before core collapse), a new 
equilibrium has been reached between the thermal pressure in the hot bubble and the ram pressure of the wind. The fragments of 
the shell dissipate into the surrounding material and the turbulent movement in the bubble starts to disappear. By the time the star 
once again reaches critical rotation, the eff'ects of the first critical rotation episode will be gone. Because the amount of material that 
has to be photo-ionized has increased during the period of critical rotation as well as because of the higher density loss rate of the 
wind, the Stromgren radius (dashed white line) has nearly receded to the edge of the shocked wind material. 
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Fig. 9. Mass loss rate and terminal wind velocity as function of time, for our 20 M© star, during the last 2 000 years before core- 
collapse. The star approaches critical rotation and the mass loss rate increases considerably. The wind velocity first decreases as 
the period of rapid rotation starts, than it increases again as the star starts to shrink. The continuous lines show the angle averaged 
values, similar to Fig.|2l The dotted lines show the mass loss rate and velocity along the polar axis, while the dashes lines give the 
values for the equatorial plane. 
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Fig. 11. Similar to Figs.[5]and|7l The final evolution of the circumstellar medium close to the star. These frames show the density of 
the circumstellar medium 21 yr and 1 yr pre-core collapse. The star is rotating at critical velocity, creating a disk of matter around 
the equator. The wind velocity has decreased, which creates a rarefaction wave that moves away from the star (R ~ 23.2.1 pc on 
the left, R ~ 3.0..3.5 pc on the right). On the inner side of the rarefaction wave, a shell is forming due to the increased mass loss 
rate. At larger radii the disk starts to expand due to its internal thermal pressure. Small changes in the escape velocity of the star 
cause a series of small waves in the wind material. The Stromgren Radius of the star lies outside the area of this simulation (Fig. [8]). 



